Schrodinger equation for the one-particle density matrix of thermal systems: An 
alternative formulation of Bose-Einstein condensation 
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We formulate a linear Schrodinger equation with the temperature-dependent potential for the one- 
particle density matrix and obtain the condensation temperature of the Bose-Einstein condensate 
from a bound-state condition for the Schrodinger equation both with and without the confining 
trap. The results are in very good agreement with those of the full statistical physics treatment. 
This is an alternative to the Bose-Einstein condensation in the standard ideal Bose gas treatment. 
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I. INTRODUCTION 

Based on the Bose-Einstein statistics^, Einstein 
preidicteda that for a gas of non-interacting, massive 
bosons, below a critical temperature a finite fraction of 
the total number of particles would occupy the lowest- 
energy single-particle state forming a condensate. This 
phenomenon of Bose-Einstein condensation (BEC) of di- 
lute gases has been extensively studied theoretically 3 
and only recently it has been experimentally^^ real- 
ized for dilute atomic gases of 87 Rb, 7 Li, 23 Na. It is 
known that the Gross-Pitaevskii equation (GPE) 7 -, a non- 
linear Schrodinger equation for the macroscopic wave 
function of the BEC provides a detailed description of its 
ground state spectrum. Though the GPE has been solved 
variationally 8 ' 9 , in this paper, we give a new theoretical 
treatment of the BEC by incorporating the symmetry of 
the system with a temperature-dependent potential and 
solving a linear temperature dependent Hamiltonian vari- 
ationally. By incorporating a non- linear term like ~ \yj\ 2 
in our Hamiltonian one can construct a finite tempera- 
ture extension of the Gross-Pitaevskii equation 7 - for an 
interacting bose gas. The present approach can be ex- 
tended in many ways including the study of temperature 
dependence of the energy of a vortex in a Bose gas and 
effect of temperature on the interference of two macro- 
scopic Bose-condensates. The free-particle many-body 
Hamiltonian^- with the quantum correction due to the 
Bose/Fermi statistics of the system is given as a function 
of temperature as 



\ - P 



(1) 



where pi and fi are the momentum and position of the 
particle, the potential is given by 



v(r) = — /cgTln 



1 ± exp[ —2ir— 



(2) 



and A is the thermal wave length, A = W '277 h 2 /mkBT '. 
The plus sign is for bosons, and the minus sign is for 
fermions. The Hamiltonian (1) represents first quantum 



correction to the classical gas high temperature expan- 
sion and describes the system in the range where the inter 
particle separation is much greater than the correspond- 
ing thermal wavelength, |rj — fj\ >> A. The statistical 
potential v(r) is attractive for bosons and repulsive for 
fermions. The one particle Schrodinger equation is given 
as 



2/y 



■V 2 V> + vip = E'tjj 



(3) 



where fi — m/2 is the reduced mass of the two particle 
system when the Hamiltonian Eq.(TT]) is written for i = 2. 
The exact thermal one-particle density matrij*ii for bose- 
systems, in position representation is given as 



p 1 Q?-?\,T) = F(\r-?\,T) 
T3 53 j ;~*exp(-aj - irr 2 /j\ 2 ) 



(4) 



where a has the value < a < 1. Our claim is the leading 
term in the one particle density matrix^i of a bose-system 
in the position representation can be approximated as 
pf(r,T) « \ip(r,T)\ 2 where L denotes the leading term. 
We replace 1-particle Hamiltonian ^-V 2 + v by 



H, 



eff 



2ft 
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ev 



(5) 



(6) 



and we have approximated v(r) « v(r) for large T, that 
is when the intcrparticlc separation is much greater than 
the corresponding thermal wavelength, r >> A. We have 
multiplied the potentail with e to take into account the 
higher order correlation effect, and it will be determined 
dynamically such that the eigenvalue of the Hamiltonian 
is < 0. We derive the dependence of critical condensation 
density on temperature of the system with or without a 
confining trap. We state an alternate derivation of the 
quantum potential and analyse our wavefunction in the 
limiting case in the Appendix. 
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FIG. 1: Variation of ground state energy of the bose gas with 
the parameter 8 for different values of e. For e = 3.89711, 
the curve with the solid line shows that the ground state just 
begins to form. The long dashed line is for e = 3.79711 and 
the small dashed line is for e = 3.99711. 



II. SCHRODINGER'S EQUATION IN 3-D 

The time independent 3-D Schrodinger equation in 
spherical polar coordinate for the reduced mass fi = m/2 
is given as 



1 



d_, sine _d_, + L_ 

2/i r 2 dr dr r 2 sin9d9 89 r 2 sin 2 
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+i(r, T) 



#^0,0) =£ty(r,0,<£)(7) 



where E is the energy eigenvalue corresponding to the 
wavefunction. We write the total wavefunction ip(r, 9, 4>), 
where T is a parameter in the Hamiltonian interpreted as 
temperature, as a product of the radial and the angular 
parts. 



iP(r,9^) = R rh i(r)Y lim (9,cj ) ) 



(8) 



Substituting the above ip(r, 9, <fi) in the 3-D Schrodinger 
equation, we get the equation for R n .i{r) as 



3^ l^ + l£l E -^T)- \.J }R(r) = (9) 



1 8 2 8R 2n 

r 2 dr dr H 2 2/j.r 2 



Using R(r) = we get the standard 1-D radial 

Schrodinger equation as 

+ [v{r, T) H — - — ]x = Ex, 



2fi dr 2 



2fir 2 



III. BEC WITHOUT HARMONIC TRAP 
A. Variational calculation 

First we study the simple case of the BEC without any 
confining trap. For I = 0, the ID Hamiltonian is given 
from Eq. (fTC)| as 



-h 2 a 2 

~2iTdr 



H, 



eff 



(10) 



We choose the trial wave function x{f) as 
X {r,T) = Are~ Sr2 



(11) 



where 5 is a variational parameter which will depend on 
temperature T and A is normalisation constant given as 



A=( 4 

7T 



(12) 



Now we calculate the expectation value of the Hamilto- 
nian H given by Eq. (|10|) , which gives us the energy 



E =< x(r,T)\H\ X (r,T) > 
We get energy E as 



E = k R T 



36 
2ixa 



1 



1 + 7rai5 



(3/2)-, 



(13) 



(14) 



where a 



mk B T 
TrrW 



and S/a is dimensionless. We note, 



even if temperature appears in a complicated way in the 
interaction term in Eq. ((6|) , the effect of T in the energy 
expression [Eq. (Tl"4")) ] is multiplicative. Now we do the 
minimization of the energy given by the Eq.(|14p with 
respect to 8, which is aE gg'^ — 0. We get 



A _ e 2/5(J_)V5 

net net 



1 = 



(15) 



Since the bound state starts to appear when the total 
energy E — > 0~ and the above minimum condition is sat- 
isfied, we therefore solve the two simultaneous equations 
for the two unknowns 6 and e as 



E(e,S) 
dEje, S) 
86 











(16) 
(17) 



and we get e = e c = 3.89711 and 5 = S c = 0.57ra. In 
the next section we will see that this set of values of e, S 
gives the right relation between density and temperature. 
The plot of E(S) as a function of S for different values 
of e is given in Fig.l. This shows for e — e c = 3.89711, 
6 = S c = 0.57ra the bound state just begins to form which 
is onset of BEC (shown as solidline). 



B. The Density- Temperature Relation 

We now calculate the density of an uniform bose gas in 
a volume v at temperature T at which the condensation 
occurs. We have effectively one particle in volume v. We 
use the standard definition of volume as v —< r 3 > ie, 

Jr 3 \ij(r)\ 2 d 3 r\ _ ( J r 3 \ X {r)\ 2 dr 



J\ip(r)\ 2 d 3 r 



J\ X (r)\ 2 dr 

= ^nP' 



(18) 
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we get 



n* = ~ = 2-46740 fe£ 



(19) 



We see from Eq. (|T9|) that we reproduce the right relation 
between condensation density n c of the bose gas and the 
temperature T at which the condensation occurs from an 
effective one particle picture. The exact relation has the 
constant^ 53/2(1) = C(3/2) = 2.612, where ({x) is the 
Riemann zeta function of x. The density-temperature 
relation is the most important result of our paper. Since 
ours is a variational one, the coefficient is little off from 
the exact value. This shows that our formalism repro- 
duces the relation which is known otherwise for an ideal 
bose gasi£. 



IV. BEC IN THE HARMONIC TRAP 

Now we study the BEC in an isotropic harmonic trap. 
The Hamiltonian for the quantum bose gas in the trap is 



given as 



H = 



-h 2 d 2 

2/x dr : 



— ekuTi 



(20) 



We have added the extra confining potential with fre- 
quency w to our original Hamiltonian [Eq. (|10[) ] and we 
want to study its effect on the BEC. We use our trial 
wave function given by Eq. ljlip with the usual normal- 
ization constant A given by Eq. (fT2")l . We calculate the 
energy E and minimize it to get one more equation to 
solve for the two parmeters e and S. The two equations 
are given as 



E = k R T 



36 
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(1 + nctS- 1 )'- 



3V2, huj ., 



(21) 



dE 
~d6 



3 e3(7ra/<5) 2 3^ Hv ^ 2 

2to 2to(1 + 7TCW5" 1 )! 8to k B T 5 

(22) 

Now we use this two equaion with E = and §f = 
to solve simultaneously for the two unknowns e and 5 at 
different values of temperature T while keeping the value 
of frequency u) constant. Since now we can not solve for 
S in a closed analytic form as a function of T, we do the 
calculation numerically. For each given value of temper- 
ature T, we solve for e and S. Then using Eq. (|18p .we 
calculate volume v and hence condensation density n c 
for each value of 5 for a set of values of T. We show in 
Fig. 2 the numerical data obtained for n c vs T and the 
fitting of the data with a function 0.354T 3 / 2 , where we 
have taken Tuv/ks = 1- The fitting is amazingly good. 
From the ideal bose-gas treatment the value of the co- 
efficient is 0.167 instead of 0.354. Since ours is a vari- 
ational one, the coefficient is overestimated by a factor 



of 2. Otherwise our formulation correctly reproduces the 
density-temperature relation for the bose gas confined in 
a harmonic trap. Since we know that siz&i 2 - of the ther- 
mal cloud is given by r ~ (7^37)^1 using this in Eq. fTO]) 
i.e, N~ v x T 3 / 2 , we get the number of particles N vary 
as T 3 . This is the relation one gets if one uses other 
methods- 2 for the calculation. This shows that the sin- 
gle particle representation of the BEC is very good. 
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FIG. 2: Variation of condensation density with temperature. 
The +'s are the numerically calculated data and the solid line 
is a fit function given as 0.354T 3 ' 72 . 



V. RELATION OF OUR TRIAL 
WAVEFUNCTION TO THE EXACT 
ONE-PARTICLE DENSITY MATRIX 

For an ideal Bose-Einstein gas it has been shownii that 
the one particle density matrix is given by Eq.(j4]). Since 
< a < 1 , we can see that the leading term in the above 
density matrix with j = I has the dependance on r and 
T as 



r) ~ ^e- 



(23) 



We can also calculate the one particle density matrix 
from our trial wave function Eq. (|lip as 

47T _x,- 2 



(24) 



Comparing Eq. (|23|) and Eq (|24]) we see that the functional 
dependence on r and T of both the density matrices are 
same. Infact Eq.([2~4f is the density matrix 13 for a three- 
dimensional system with N particles. So it justifies the 
one particle description of the BEC in our formalism by 
taking a temperature dependent pure state. 



VI. CONCLUSION 

In conclusion, we have shown that it is possible to 
get Bose-Einstein condensation with or without the har- 
monic trap by considering a Hamiltonian with symmetry 



4 



dependent potential and a temperature dependent pure 
state. We reproduce the right relation between condensa- 
tion density and temperature at the critical point. This 
is an alternative way to the Bosc-Einstein condensation 
in standard ideal bose gas treatment. By incorporating 
a non- linear term like ~ \ip\ 2 in the linear Schrodinger 
equation one can construct a finite temperature exten- 
sion of the Gross-Pitaevskii equation^. The present ap- 
proach can be applied to study the temperature depen- 
dence of the energy of a vortex in a Bose gas and effect 
of temperature on the interference of two macroscopic 
Bose-condensates. 



APPENDIX A: DERIVATION OF THE 
INTERACTION POTENTIAL FROM THE 
DISTRIBUTION FUNCTION 

For a pair of non-interacting particles in a volume v 
the canonical partition function is given as 



(Al) 



where A is the thermal de-Broglie wave length defined be- 
fore. It can be shown- - that the probability distribution 
for the separation r is given as 



which is a kind of entropy which is given as 



S = — kn In 



1 ± exp(-2nr 2 / A 2 ) 



(A4) 



Now if we multiply the entropy S with temperature T 
we have energy, which here appears as the interaction 
potentail energy. 



V = ST =-k B T\n 



I ± exp{-2irr 2 / A 2 ) 



(A5) 



This is exactly the interaction potential^ used in our 
theory. 



APPENDIX B: BEHAVIOUR OF THE WAVE 
FUNCTION WHEN T -> oo AND r -» 

The interaction potential Eq. Q in the limit of high T 
behaves 1 ^ as 



lim v(r,T) = —T 4 5(r) 



(Bl) 



And our trial wave function Eq. lfTTj) behaves^ 5 - in the 
limit of high T as 



Pl2 



1 



1 ± exp(-2irr 2 / A 2 ) 



(A2) 



Now we take the logarithm of the above equation and get 



In Pi 



12 



hi 



1 ± exp{-2irr 2 / X 2 ) 



(A3) 



where +(-) sign is for boson(fermion). In writting the 
above equation we have dropeed the constant term in- 
dependent of r. We multiply by — kg to find a quantity 



lim x(r, P) = cT4-r5(r) 



(B2) 



where c is a function of h and m and independent of 
r and T. If we use the high T limit of interaction po- 
tential Eq.(Bl) in the 1-D radial Schrodinger Eq. (fT0|) . 
the above hight T limit of \( r :T) [Eq.(B2)] does satisfy 
the equation with an eigen energy which is not bounded 
from below. The eigen value diverges. This is exactly the 
behaviour of the variational energy valuc[Eq. (|14[) ] which 
also diverges in the high limit of T. 
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